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Abstract. A deterministic coalescing dynamics with constant rate for a particle system in a finite 
volume with a fixed initial number of particles is considered. It is shown that, in the thermodynamic 
limit, with the constraint of fixed density, the corresponding coagulation equation is recovered and 
global in time propagation of chaos holds. 
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1 Introduction 

Coagulation equations are widespread models describing large sets of coalescing particles. It is argued 
in such descriptions that the coagulation rate depends on the masses of the particles. In the spatially 
homogeneous models it is further assumed that spatial fluctuations in the mass density are negligible. 
In the first example of this sort, proposed by Smoluchowski [21], particles of radius r moving by 
Brownian motion with variance proportional to 1/r meet at rate proportional to (ri +r2)(r^^ +r^^). 
The concentration f{t, m) of particles of mass m = 1, 2, 3, • • • at time t satisfies the equation: 

-^{t, m) = -Y, K{m - k, k)f{t, m - k)f{t, k) - f{t, m) J] K{m, k)f{t, k) 

k=l k=l 

K{m, k) = (m^/^ + (m~^/^ + k~^/^^ 

If masses are considered to take continuous positive values we are generally led to the continuous 
version: 

—-{t,m) = — K{m — m' ,rn')f{t,m — rn')f{t,m')drn' — f(t,rn) K{m,m')f{t,m')dm' (1-1) 
2 Jq Jq 

Different coagulation kernels K{m, m') may be considered, in particular constant, additive and 
multiplicative kernels have been widely studied. The concentration /(t, m) is defined as the average 
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number of clusters of mass m per unit volume at time t in the discrete case, and as the average number 
of clusters of mass in [m,m + dm] per unit volume at time t in the continuous case. Our purpose is to 
deduce the continuous coagulation equation as the limit of equations describing a coalescing system 
of finite (large) number of particles in a finite spatial volume by means of PDE arguments. 

To this end we consider, in a given volume V, a system of particles whose number changes in time 
due to the coagulation process. At any time t > 0, given any € N* and (mi, m2, • • • ,miq) € (M^)^ 
we consider the mass distribution functions PAr(t, mi, ttitv)- Each of these functions describe the 
state of a system constituted by N particles in a volume V where, at time t, one particle has mass 
between mi and mi + dmi, one has mass between m2 and m2 + dm2 and so on. We may then define 
the probability to have N particles in the volume V at time t = as the following: 

/'OO /'OO 

P{{),N) = — I ■■■I P%{mi,...,mN)dmi---dmM 
■ Jo Jo 

and 

oo 
N=l 

Given any Nq G N* we consider initial data satisfying the following conditions: 



' (i) P^(mi,...,?nAr) = 0, for ah A /Ao, 

(ii) P^^{mi,...,mNo) = {Noy.fo{mi) (S> ■ ■ ■ <Si foiruNo), mi€(0,+oo), i = 1,2,. . . Nq ^-^ 
/•oo 

where: /o > 0, / fo{m) dm = 1. 
Jo 

Condition [i) expresses that at t = the system has exactly Aq particles. It is easily seen that 
conditions (i) and [ii) imply P(0, A) = 5{N — Nq). The average number of particles is then: 

oo 

^ AP(Ar,0) = Ao. 

Af=l 

as expected. 

The starting point of our analysis are the evolution equations satisfied by the mass distribution 
functions P/v(t,mi, ...,m]\j) throughout the coagulation process. These equations had been obtained 
in [IS] and later in [H], (22] and [H]. It is then possible to deduce a system of equations for the 
correlation functions: 

fj{mi, . . . ,mj,t) = yZjM ^ dmj+i.../ dmM PN{mi, . . . ,mN ,t). 

jSf^j v-'^ ~ J>- Jo Jo 
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where j € N*. Notice that for j = 1 we have 

Pn {t, mi , 1712 ) • ■ ■ ,^n) dull • • • dm n 



fi{t,mi)dmi 

and therefore fi{t) is the density function associated to the mean number of particles at time t. 
Consider now the fohowing rescaled functions: 



fj {t,mi,--- ,mj) = ^ — (1.3) 



We have, for j = 1: 



/ fY i^i^i)dmi = — — — = mean number of particles per unit volume at time t>0. 
Jo ^ 

We then study the limit of these rescaled functions as 

A^o ^ +00, V +00, -y po. (1.4) 

By (ILII), we have 

//(O, mi,..,m,) = mNo-l).^iNo-J + l) ^^,^^^^ . . . ^ 

and it easily follows that: 



Vj€N*: hm Wf^ (0) - (poY f^'lk^iR+y) = 0- (1-6) 

Nr. 



The main result of the paper is the following. 

Theorem 1.1 Suppose that the coagulation kernel K is constant. Let {/|^(t)}jgN* be the sequence of 
functions defined in il.3\) with initial data il.5\) . Then, 

Vj>l: hm ||//(t)-/(t)®^||Li((R+).)=0 (1.7) 

Nr. 

where f is the unique solution in C([0, +oo); L^(M"'")) of the coagulation equation il.l]) with initial 
datum po/o- 

The convergence result ()1.7p under the hypothesis (jl.Sp for the initial data is usually known as propaga- 
tion of chaos. By the initial condition (jl.Sp the particles are identically and independently distributed 
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at time t = 0. The conclusion ()1.7p means that at any time t > 0, the system of particles has still 
that property but only asymptotically in the limit ()1.4p (cf. [18] and references therein). 

The suitable way to describe the finite system of coalescing particles is the stochastic Markov 
process called Marcus-Lushnikov process. This has been introduced in [15] and later in [8], [22] and [14j . 
In these four references the authors obtain and pay special attention to the evolution equation satisfied 
by the probability of the stochastic process, called joint frequency or mass distribution function. They 
deduce from that equation the system of equations satisfied by the moments of the mass distribution 
function. The conditions under which this system may be approximated by the coagulation equation 
are discussed. It is seen in particular that this depends on the no correlation between the numbers of 
particles of different mass (cf. also |20j for the continuous case). 

In [To], [llj . [17j . [7] the authors consider the Marcus-Lushnikov process itself. The weak con- 
vergence of suitably rescaled versions of these processes, called stochastic coalescents, towards the 
solutions of the coagulation equation is proved under different conditions on the coagulation kernel 
K, both in the discrete and the continuous case. 

We believe nevertheless that the questions raised in [15], [8], [22] and [H] about the system of 
equations satisfied by the moments are of interest. In particular, to study the case of continuous masses 
via the equations for the mass distribution functions by means of PDE arguments is worthwhile, for 
itself and for further work [B]. That leads very naturally to a set of equations for the correlation 
functions /j's that, as we have already said, has been discussed to some extent in [15], [8], [22] and 
[14j . The system of rescaled correlation functions is similar to the BBGKY hierarchy appearing in 
the study of many particles hamiltonian systems. Notice that in such cases the underlying microscopic 
dynamics is deterministic (since it is given by the Newton equations associated with the Hamiltonian 
under consideration), while the dynamics of the coalescing particles is probabilistic, given by the 
Marcus-Lushnikov process. Nevertheless, our system of countably coupled equations may be treated 
along similar lines and our result may be stated in the same language in terms of propagation of 
chaos. These arguments inspired from many particles hamiltonian systems have already been used 
in the context of particle's coalescence in |13j where propagation of chaos is proved for a coagulation 
equation with spatial dependence and constant kernel. 

Let us recall here that for mean field models of classical particles with suitable two body potential, 
the propagation of chaos holds true [1], [3, [l6] and may be seen as a law of large numbers. The method 
of the proof is based on the use of the so called empirical distribution and, as explained in [S], it relies 
upon the two following facts: the empirical distribution is a weak solution of the Vlasov equation; 
weak solutions to the Vlasov equation are continuous with respect to the initial datum in the topology 
of the weak convergence of the measures. In that case no proof of this result based on PDE arguments 
is known (cf. for example the discussion in Section 1.4 of |18|). 

When the coalescing kernel is constant, as it is considered in this work, the proof of the main 
result is rather simple and clear, due to the fact that the suitable functional space, the set of non 
negative integrable functions over M"*" , is easily seen to be globally preserved along the time evolution of 
the system. The stochastic coalescent for K = 1 is known as Kingman's coalescent. The construction 
goes back to Kingman [T2] and it has been extensively studied. Some details and references may be 
found in [2] and [3] . For more general kernels the proofs may become more technical and involved and 
some other delicate questions may appear (cf. [B]). 

The plan of the paper is the following. In Section [2] we briefly recall the description of the finite 
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particle system and present an existence and uniqueness result for the system of equations satisfied 
by the mass distribution functions. In Section [3] we introduce the correlation functions and deduce 
the system of equations that they satisfy. Two simple time asymptotic properties of that system are 
derived in Section HI The rescaled system is introduced in Section [5] where the main theorem is proved. 



2 The particle system 

We want to describe a finite system of particles contained in a finite volume V whose number N is 
not fixed through time. To this end we consider the probability space defined as follows. 
The sample space is: 

oo 

= (J (iV, C N* X ^(M+) (2.1) 

N=l 

^Ar(M+) : subsets of R+ with cardinal N. (2.2) 

It is the set of all the pairs {N,uj]y) where iV € N* and a;jv is any finite subset of N positive real numbers 
"T'l; "1-2 5 ■ ■ ■ 5 miyf. For every time t > the probability distribution of each state {N, (mi, . . . , mjy)) € 
is Pjyit, mi, . . . ,mN)/N\ where {PN{i)} Nizfi* is a sequence of non negative functions PN{t) = -P/v(*5 ■)) 
each of them defined on and normalized according to: 



CO ^ 

N=l 



/•oo /'OO 

/ dmi ■ ■ ■ (imArP/v(t, TTT-i, • • • , mAr) = 1 (2.3) 
Jo Jo 



The functions P]\fit,mi, . . . ,mAr) are assumed to be symmetric with respect to any permutation of 
the indices 1, . . . , and no restrictions are imposed on the range of mass values other than > 
for i = 1, . . . ,N. Each P/v = -P/v(i)'^i) ■ ■ ■ tItt-n) is the mass distribution function of the A^-particle 
configuration (mi, . . . ,m]\f). For every € N* the function: 



CO 



^'(t,A^) = ^y dmi---J dmNPNit,mi,--- ,mN) (2.4) 

is the probability that at time t the system is constituted by N particles and the normalization ()2.3p 
is natural. 



Remark 2.1 The factor 1/N\ in definition \2.4^ is needed to compensate for counting all the N\ 
physically equivalent ways of arranging the masses mi, . . . ,m7v in order of size. 

Let ^y(m, /i) > be the rate (per unit time) at which a single particle of mass m collides and 
coalesces with a particle of mass fi (given only that each is in the considered spatial region of volume 
V). We assume that: 

Av{rn, jj.) = AvifJ-^rn), (2.5) 
as it is reasonable from a physical point of view. 
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In a coalescing system the number of particles will be varying along the time evolution (actually 
it will be decreasing) so what is meaningful is to consider the time evolution of the sequence of mass 
distribution functions {P]\f{t)}]y^^* . As explained in [15], [8], [22] and [H] the evolution equations for 
the mass distribution functions are: 

dtPN = dfi Avirui- fi,fi)PN+i{mi,...,rne-i,me- fi,me+i,...,mN,f^,t) + 

1 ^ 

- -^^PNimi,. . . ,mN-,t)Av{mi,mq). iV=l,2, ... (2.6) 

The first term on the right hand side of ()2.6p is the so-called gain term and it describes the positive 
contribution due to the coagulation of a particle of mass fi, with n < m^, with a particle of mass 
rrii — (in an + 1-particle configuration), giving rise to a particle of mass mi (in an A^-particle 
configuration). The second term on the right hand side of (|2.6p is the so-called loss term and it 
describes the negative contribution due to the coagulation of a particle of mass mg with a particle of 
mass niq (in an iV-particle configuration), giving rise to a particle of mass nig + niq (in an — 1- 
particle configuration). That set of equations completely neglects the contributions that would be due 
either to many-boby collisions (e.g., three-body collisions passing from an A^-|-2 configuration to an A^ 
configuration or from an A^ configuration to an A^ — 2 configuration ) , or to the occurrence of multiple 
binary collision (e.g., double binary collision again passing from an A^ + 2 (or A^) configuration to an 
A^ (or N — 2) configuration). 

An important feature of the system (|2.6p is that the total mass is preserved. More precisely, de- 
noting by Mat the total mass of the configuration mi, . . . , m^r, i.e. mi + ■ ■ ■ + m]\f = Mjy, then the only 
processes that have a non zero contribution in the time variation of the distribution Pi\f{mi, . . . , rriTv, t) 
are those associated with configurations with the same total mass Mtv- In fact, the gain term in ()2.6p 
takes into account A^ -|- 1-particle configurations with total mass mi -|- • • • -|- mg — ^ + ■ ■ ■ + m^ + A* = 
Mat and the loss term in (j2.6p leads to A^ — 1 particle configurations in which the total mass is 
fni-\ h {mt + mq) H h m^ = Mn- 

The system of equations ()2.6p has been studied in [15], [8], [22] and |14j . For the sake of com- 
pleteness we present in this Section an existence and uniqueness result that suits our purposes. 



2.1 The coalescence rate 

In general the coagulation kernel Av{mi^m2) at which two particles of masses mi and m2 coalesce in 
a volume V , has the following form: 

Av{'mi,m2) = Cv A{mi,m2). (2.7) 

The function A{mi,m2) encodes only the dependence of the coagulation rate on the masses mi and 
m2- The term Cy contains the dependence of that rate with respect to the volume V as well as that 
on other physical properties of the coalescence process under consideration. An example (considered 
for instance in [15] and [8]) is the following: 

Cv = ^E{\vi-V2\)7r (2.8) 
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where E{\vi — V2\) is the average relative velocity of the two particles. The volume dependence of 
Cy takes into account the fact that the coagulation rate increases as the proportion of the volume 
occupied by the particles with respect to the total volume V increases. That gives a dependence 
inversely proportional to V, i.e. like 1/V. We consider in this paper the simplest possible case 
compatible with our purpose, that is to study the limit of the finite particle system as ^ ^ +00 at 
fixed density, namely: 



Av{mi,m2) = ^ (2.9) 



Other kernels are considered in [6]. 



2.2 The initial data 

We consider system (j2.6|) with initial datum: 

r (iVo)! Mmi) . . . MruNo) , if iV = iVo, 
P^(m^) = <^ (2.10) 
[ if iV / iVo, 

where the function /o is such that: 

/•oo 

/o(m)>0 a.e, / o!m/o(m) = 1. (2.11) 
^0 

By choosing the initial datum as in ()2.10p - ()2.1ip . and defining: 

P'^{N) = — I dmi... I dniN P%{mi,...,mN), (2.12) 
for every € N*, we have 
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P^{N) = 1 (2.13) 

N=l 

and condition ()2.3p is satisfied at time t = 0. With that choice of initial datum: 

P°(iV) = (5(iV-iVo), 
i.e. at t = our system has exactly Nq particles. 

2.3 Well-posedeness of the equation for P/v(t) for the constant kernel 

For every G N* we define the operator Gn, mapping N + 1-particle functions into A^-particle 
functions, as follows: 
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GN[QN+i]{mN) = -^/^ duiN+i QN+i{mi 

£=1 



N 



(2.14) 

where Qn+i € L^((]R+)^"*'^) and we denote: 

uiN := (mi, . . . ,mN). (2.15) 
We shall also use the following notation: 

/ft ftN-l 

dtN ■■= / dh... dtN- (2.16) 
Jo Jo 

and denote Ti the set: 

-H = {{fNjNm^; fk G i'((K+)'), k G N*} (2.17) 

Definition 2.2 We say that a sequence {-PAr(i)}AreN* S V. solves the system 112. 6\) if, for every t > 
and every N G N*, each term in ( fO) belongs to L^{{R+)^) and the equality holds in L^{{R+)^). 

Lemma 2.3 For every Qn+i G L1((M+)^+1); 

N 

\\Gn[Qn+i]\\l^(^(^^+)N) < — ||(57V+l||j;^i((]8+)iV+i) ■ (2.18) 

If Qn+i G L'^{{M+)^^^) and Qn+i > then 

N 

\\Gn[Qn+i\\\l'L({r+)i^) = 2Y II^^+iIIli((k+)^+^) ■ (2-19) 
Proof of Lemma 12.31 We start proving ()2.18|) . To this end we write: 



Gn[Qn+i]\\li((r+)N) = dmi--- diUN \GN[QN+i]{mi,- ■ ■ ,mN)\ 



-j^ ^ /"OO /"CO 



diUN dmN+i\QN+i{mi,m2,--- ,me-mN+ir-- ,mN+i)\ (2.20) 



•me 

dmi I dmN+i\QN+i{'mi,m2, - ■ ■ , - mAr+i, • • • ,mN+i] 



where in the last step we have used that the integrals with respect to the variables mj for i = 1, - ■ ■ ,N 
are over independent domains and therefore can be performed in any order. The final and trivial step 
is to use Fubini's theorem in the integrals with respect to and mN+i- That may be done since by 
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hypothesis Qn+i G L"^((M'^)^'^"^) and therefore, for ahiiost every {mi,m2, ■ ■ ■ , m£_i, m^+i, • • • ,m]\f) 
QN+i{mi,m2, ■■■ ,mi,--- ,mN,mN+i) G L^((R+)2). 

If we also had Qn+i ^ the inequahty in (|2.20p would be an equality and we would have ()2.19p . 

□ 

We can state now the existence and uniqueness result of solutions for system ()2.6p . 



Theorem 2.4 Suppose that the Ay = V""^. Then, the Cauchy problem for the system h2. b]) with 
initial data {P^jjvgN* defined in /12.10\ )- [2.11\) has a solution {-P/v(t)}7VeN* £ ^ for any t and such 
that Pn € C°°([0, +00), L^((]R+)^)) for every N € N*. Such a solution is given by 



PNit,mN) 



j dtNo-N {e-^^ Gn}... 



(iV0-l)(iV0-2) 



e *'Vo-ivpO^(m^j, /or 1 < iV < iVo 



(iVo)(iVo-l) 
e 2V- 



*pO,^(m^J, forN = No 

[ 0, for N > iVo. 
Moreover, for every t > 0, -PAr(t) > and: 



(2.21) 



with 



N=l 

P(t,N) := ^ 
^ ' ^ AT! 



dmi 



dm]\fPj\f{t,mi, ■ ■ ■ ,miy) 



(2.22) 
(2.23) 



That solution is unique in the set of sequences {PN{t)}N^N* G 71 such that P n G C^([0, +(X)); L^((R"'")^)), 
Pat > for all N and the function hp{t) := E^=i ^i^, N) satisfies hp G Ci([0, +00)). 

Proof of Theorem [23] By our hypothesis P^^ G L^((M+)^o). By Lemma ([2:3]) we deduce that 
for every 1 < < Nq: 

Then, for any a and r, the function 



satisfies: 



r,{t, mjv) = e'^(*-^)GAr Gn+i ■ ■ ■ G^o_i[P^J(mjv) 



^ec-([o,+oo), lH(m+)^)). 



It then follows that PAr(t) given by ((2:211) is such that Pn G C°°([0, +00), L1((M+)^)). 

We check now that the above sequence {PN{t)}Nm* satisfies the system (j2.6p . The equations 
for N > Nq + 1 are trivially satisfied. Suppose now that N < Nq — 1. Taking the time derivative of 
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the left hand side of ()2.2ip we obtain: 

dPN, , N{N-1)^ , 
-^(t,mjv) = — P7v(t,mjv) 



dt2 / dts - . . 

^0 ^0 



dtNo-N Gn {e-^^ (*-*2)G;vo-i} • • • 



N(N-l) 



G 



N 



2V 

/■* Z"*! ftNo-N-2 f (N+1)N . . 

dti dt2... dtNo~N~l [e — ^'~''>Gno~i\ ■ ■ ■ 



And, since 



* /■*! ^■tNo-N-2 ^ f (iV+l)iV 

dti / dto . ■ 



we obtain 

at 2v 
that is nothing but the A^^'^ equation of system (|2.6p . If A = Aq, we have: 
dPNo., . Ao(Ao - 1) „ 



9t ' ' 2V 

NojNo - 1) 

2y 



PAr(t,mAro) + GnoIPno + iWi'^No) 



since, by ()2.2ip . P/VQ+i(t) = for any t. 

We check now ([2^. To this end we first notice that, by ([2:21]) : 

drriN Pn (t, mi , . . . , mAr) . 



N=l N=l •^'^ 
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Therefore, using ([22]) and ([2T9]) (since Pn > 0): 

d °° N( N — 1) ^0 

-j;P(i,iV) = -E N^2V ll-P^(^)II^M(M+r) + E]^ll^^[-P^+i]||Li((M+r) 

Af=l Af=l ■ Af=l 

^iV(iV-l),, N ,, 

Ar=l ■ Af=l 

^iV(iV-l),, ,,,, ^ M-1 ,, ,,,, 

Af=l ' M=2 ^ 

^iV(iV-l) ^M(M-l),, , ,,, 

N=l ' M=l 
Identity (|2.22|) follows from the fact that, by continuity and (j2.13|) one has: 

No No 

^Pit,N) = ^P\N) = l. 

N=l N=l 

In order to prove the uniqueness, we first show that any solution {PNi't)}NeN* , such that PN{i) > 
and hp{t) = Y1n=i ^(*' ^) verifies h-p € C^((0, +oo)), satisfies Pjv(t) = for ah t and iV > A^o + 1- 
To this end we write: 

d °° °° N( N — 1) °° 

- E ^(*'^) = - E N12V E \\GN[PN+l]\\L^iiR+)N) 

N=No+l N=No+l ' N=No+l 

^ N{N -1) - ^ N - 

- ~ 1^ N\2V ll^^^*^llLi((R+)iV) + 2^ ^y^ll[^^+lJllLi((R+)iv+i) 

N=No+l ' N=No+l 

^ Ar(iV-i) ^ iV(Ar_i) 

= " 2^ Ari9F + 2^ ^,2y ll[^ivJ||Li((R+)iV) 

Ar=Aro+l N=No+2 



Yy|^ll^^o+lWllLi((K+)^o+i) < 



We deduce that: 



E P{t^N)< E nO,A^)= E ^°(^)- (2.24) 

Ar=Aro+l Ar=Aro+i N=No+l 

Since, by hypothesis, the right hand side of ()2.24p is zero and PN{t) > we deduce that -PAr(t) = 
for every N > Nq + 1 and t > 0. That leaves only a finite number of equations in the system (|2.6p . 
for = 1, • • • , Nq. Moreover the equation for N = Nq yields: 

„ iVo(iVo-l) 

Pno {t, mjvo ) = Pno (mjvo ) e * • 

Then, the finite system decouples and may be explicitly solved for N = I,-- - ,Nq — 1 to obtain 
expression (j2.2ip . □ 
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Remark 2.5 One would expect that for every initial data {P^^New ^ ^ such that > for 
every N and '^P^{N) = 1 there exists a solution {Piy{t)}iy^fq* such that, for every N: PN{t) > 0, 
Pn G C([0,+oo);L1((M+)^)) nCi((0,+oo);Li((]R+)^)), hp G Ci(0,+oo) and E?^=i^(*>^) = 1- 
Since such a result is not necessary for our main purpose, the proof of Theorem we only prove 



the simpler result in Theorem\2.4 



3 Correlation Functions 

For any fixed j G N* , we define the j-particle correlation function fj{mi, . . . , mj,t) at time t as 

fjjmi, . . . ,mj,t) = I dmjj^i... dniN PN{mi, . . . ,mN,t). (3.1) 
j V J ) ' J J Q 

Some well known general properties of such functions fj are the following. At any time t the expected 
(or mean) number of particles N(t) defined as 

oo 

N{t) = ^NP{t,N) 

N=l 

satisfies: 

oo 



N{t)= / fi{mi,t)dmi. (3.2) 







The function fi is then the density function associated to the average number of particles. More 
generally one may define: 

oo 

{N{t)){N{t) - 1) . . . {N{t) -j + l) = Y^N{N-l)...{N-j + l)P{t, N) 

N=l 

and one has 



(iV(t))(iV(t) - 1) . . . (Nit) -j + l) = (3.3) 

The functions fj 's will be called correlation functions since their definition is very similar to that of the 
classical correlation functions in statistical mechanics (cf. [TO]) and satisfy properties (j3.2|) and (|3.3|) . 
Notice nevertheless that they slightly differ from the correlation functions in statistical mechanics 
since in particular, there is no activity parameter (cf. |19j). 

As an immediate Corollary of Theorem 12.41 we have the following. 

Corollary 3.1 Consider the solution {P/v(t)}jvgp^. of the Cauchy problem 112. 6)) - ![27W\) - 112. 1 i]) . whose 
existence and uniqueness have been proved in Theorem \2.4\ Then, for every j £ N* the functions fj{t) 
defined as 

fjjmi ■■■ ,mj,t) = / dmj+i.../ druN PN{t,mi, . . . ,mN) (3.4) 
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are such that: 

■^^•W= /i^^^'-M fo'^ Vjg{l,...,iVo}; /,(0)=0, Vi>iVo + l. (3.5) 

(n) > V J > 1 and V t > 0; = V t > i/ j > iVo + 1- (3.6) 

(m) gC-([0,+oo);L1((IR+)^')). (3.7) 

(^^) ll/.(t)llLH(M+M)<2^~'^o! 77 E ^^^""^^"^ /«^J' = 1'--- '^0- (3.8) 

Proof of Corollary 13. II Only (|3.8|) needs an explanation. It follows from the fact that, by definition: 

No ^ 

where we used that P/v(t) = for N > Nq + 1. On the other hand, by Theorem 12.41 for eyery N < Nq 

dtN,-N / dxxvNe —^'-'^^GN...GNo-ie ^^^*^o-ivp^^(mjvo,0). 

Using ()2.19p in Lemma 12.31 we deduce: 

fNo-N 

ll^7v(t)||LH(M+n < F-(^"-^)iV(iV + l)...(iVo-l)||PM,(0)||^i((j,+).o) = 

(A^o-Af)! (A^-1)! ll^^o(U)||ii((M+)iVo)- l-^-lU) 



Using the inequality 



we obtain from (|3.1U|) : 



ll^iv(t)||Li((R+n < (AT^-AT)! ^"^"^ y-(^O-iV) 2^^0-1 (^^_^), ||P^„(0)||^,((^+).„) 

where we used assumptions (j2.10p - (j2.1ip on the initial datum P/Vo(0). Estimate (j3.8p follows from 
(I32D and (l3T2]) . □ 



3.1 The system of equations for the functions fj 

An immediate consequence of the aboye calculation and Corollary ()3.ip is the following: 
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Theorem 3.2 Suppose that {PN{t)}NeN* is the unique solution of 112. 6\) with initial data {P^}NeN* 
defined in 112. 10]) - [2.11]) and coagulation kernel Ay = Y~^. Then the set of functions {/j(i)}jeN* 
defined in solves the following problem: 

(3.13) 



dtfj = 2y'12 dn fj+i{mi,...,mi- n,...,mj,ij.,t) + 



"''^■^ fj(.TTT-i,---,mj,t)-^ I dfj,fj+i{mi,...,mj,n,t), j G N*, 



2V 2 ' 7' ^ y 

Moreover the sequence {fj{t)}jQ^* is the unique solution in C([0, +oo); L'^{{R+y))nC^{{0, +oo); 
of /]3.13]) with initial data given /]3.5]) . 

Proof of Theorem 13.21 A straightforward computation from (|2.6p using the definition (|3.4p shows 
that for any j < Nq the function fj satisfies the equation : 

No -, 



mi 

druN+i 



PiV+i(mi, ...,me- ruN+i, ■ ■ ■ ,mN+i,t) + 



1 N(N — 1) f 

~~ X] f/v - 'V W dniN,jPNimi,...,mN,t), (3.14) 

N=j ^ 

where from now on we use the notation: 

/fOQ fOO 
dm]\[j := / drrij^i . . . dmjq. (3.15) 
Jo 

The first term in the right hand side of ()3.14p . or gain term, gives the following contributions: 
^^^UTt / dmMj / dniM+i PN+i{mi,. . . ,me,- niM+i,. . . ,mN+i,t) 

(■mi 



X rmi 

= l^zZ j dmN+ifj+i{mi,...,mi-mN+i,...,mj,mN+i,t) + 
£=1 

1 ^ N-j f 



N=j 

PAr+i(mj,mj+i - rriN+i, ■ ■ . ,mN+i,t), (3.16) 

where we have split the sum with respect to i in two parts, I < i < j and j + 1 < i < N, and used the 
symmetry of Pjy with respect to any permutation of the indeces (that is preserved by the dynamics), 
to write: 

y] / druN+i PN+i{mi,...,mi-mN+i,...,mN+i,t), 

f-mj+i 

= {N - j) / dniN+i PN+i{mi,...,mj+i-mN+i,---,mN+i,t). 
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Using Fubini's theorem and definition ()3.4p we obtain 

No 



2\7 ^ (iV3j)T J J dmj+i J dmN+iPN+i{mj,mj+i - rriN+i, . . . ,mN+i,t) 

poo roo 

= TprP dniN+i / (imj+i /,+2(mi, . . . ,mj+i,mAr+i,f). (3.17) 
Jo Jo 

On the other hand, by simple algebraic manipulations, the second term in the right hand side of (j3.14p . 
or loss term, gives: 

^ 1 N(N -I) , , . 

2^ (jV _ j)! 2V io J+i-'-yQ dmNPN[mi,. . . ,mN,t) 

j(j - 1) j /■°° 

= — ^ — fj{mi,...,mj,t) + — dmj+ifj+i{mi,...,mj+i,t) + 

-j^ roo poo 

+ — y druj+i J dmj+2fj+2{mi,...,mj+2,t), (3.18) 

By (IXTiD . (IXTHD . (l3J71l and (l3J8]l . (1X13]) follows. 

Only uniqueness remains to be proved. To this end we notice that since i-V(i) = for every 
> A^O; the same is true for the sequence {/j(i)}jeN* > namely, fj{t) = for j > Nq. Therefore the 
system (|3.13|) only contains a finite number of non trivial equations and, since the equation for j = Nq 
only involves fNoi't)^ the system decouples. We conclude using the same argument used in Theorem 
12.41 to prove uniqueness for the mass distribution functions {P/v(t)}ArgN* • Q 



Remark 3.3 If the initial data /q satisfies not only f2.11\) but also: 

niQ = / mfQ{m)dm < +oo 



then, since, as we have seen above, fi{t) is the density function associated to the average number of 
particles at time t, the quantity: 



M{t)= / mfi{m,t)dm (3.19) 
Jo 

represents the average mass in the system at time t. If we integrate the equation for j = 1 in system 
a simple application of Fubini's theorem yields 



d_ 

dt 



mfi{m,t)dm = (3.20) 



and then, using \2.11\) and h3. 5|) 



M{t) = M(0) = mo A^o Vt > 0. (3.21) 

This mass conservation property is a well known feature of the coagulation equation with constant 
kernel. 
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4 The solution of the finite system for 

In this Section we study the behavior of the sequence as t — +00. 

Proposition 4.1 Let {fj{t)}^^i be the solution of ^3.13]) with initial datum ^3. 5|) given by Corollary 
CO Then 



hm N{t) = 1 and lim var{N){t) = 0, 

t— >+oo ^+00 



where var{N)[t) = {N{t))'^ — N{t) is the variance of the distribution on particle numbers P{t,N). 
The proof of Proposition 14.11 follows from the three following auxiliary results. 
Lemma 4.2 For all t > and j: 

^ll/i(OllLi((K+)^) = ~^ll/i+l(OllLi{(M+)J+i) - "^'^2^ \ \fj{'t)\\L^{{R+y) < 0. 

Proof of Lemma 14.21 Since, by definition, fj{t) > for any j and t, (the) Lemma 14.21 follows by a 
simple integration of (|3.13|) with respect to mi . . . mj . □ 



d 



From (j3.2p and Lemma 14.21 we deduce 
Corollary 4.3 (i) For all t > 0; 



4t ^ ° 



('iij There exists N^o G [1, -/Vq] such that N{t) — > A^oo as t ^ +00. 

Proof of Corollary 14.31 Only the fact N^o > 1 needs perhaps to be explained. From (i) we deduce 
the existence of A^oo G [0, A^'o], the unique limit point of N{t) as t ^ +00. On the other hand, 

N{t) = NPit,N) > Pit,N) = 1 

Af=l Af=l 

from where N^o > 1- □ 
We also deduce: 

Corollary 4.4 For all t > and every j > 1: 



(iV(t))(iV(t)-l)...(iV(t)-i + l) <e"— * N^, (4.1) 
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Proof of Corollary 14.41 By Lemma 14.21 



j(i - 1 



< 

2V 

that implies, by (|3.5p 



II/.WIIli((r+).), (4.2) 



II/.WIIlM{m+M) < e-"^* |l/i(0)ILi((M+).) <e-'^*iVi]', (4.3) 



for any time t. Since, by definition, {N{t)){N{t) - 1) . . . (A^(t) - j + 1) = ||/j(t)||ii( 
this concludes the proof □ 

Proof of Proposition 14. 1\ By Corollary 14.41 for j = 2: 



ll/2(t)llLi((R+)2) = (A^(t))(A^(t) - 1) < e~v ^0, as t ^ +00 (4.4) 

and then 



lim (iV(t))2 = lim N{t) = Noo- (4.5) 
Therefore, the yariance of the distribution on particle numbers for large times is giyen by: 



lim var{N){t)= lim {{N{t)y - N{t) ) = N^o - N^^. (4.6) 

t— >+oo i— >+oo 

Since, by definition, the yariance is non-negative 

A^oo - iV^ > ^ iVoo(l - iVoo) > ^ < iVoo < 1. (4.7) 
On the other hand by Corollary 14.31 we know that N^o > 1, thus it has to be A'^^o = 1 and therefore 

lim var{N){t) = 0. (4.8) 

t—^+oc 

□ 

Remark 4.5 In other words, the distribution on particle numbers is going to be a delta distribution, 
as it was at time t = 0, but now centered at Noo = 1- As expected, due to the coalescence dynamics, 
the system is going to be constituted by only one large particle of mass Nq ttiq . 

5 BBGKY hierarchy of the rescaled correlation functions 

Our purpose in this Section is to consider the limit where the yolume V and the initial number of 
particles A'^o go to infinity in such a way that 

iVn 

vJ™. -rr=Po€(0,+oo). (5.1) 

V, Nq^+oo V 
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To this end let us define the rescaled correlation functions as 

//(mi, . . . , m,,t) := M^^^ -^y^^^') ^ j = l,...,N,. (5.2) 

Since, as we have seen in Section [3l the function fi is the number density function, the rescaled 
function fY is the density function associated to the concentration of particles (number of particles 
per unit volume). 

By p.l3p the functions satisfy the following set of Nq equations: 

j 



= \Y. d^,fy_,,{m„...,m,-^^,...,m„^,,t) + 



- '^^^2V^^ fj ijn^i^ - ■ ■ i'mj,t) - 3 dfi fj'^^imi, . . . ,mj, n,t), (5.3) 

for j = 1,2, ... , Nq. We will refer to the family of equations (|5.3p as BBGKY hierarchy by analogy 
with the system arising in the framework of many particles hamiltonian systems. 
By (|3.5p and ()5.2p the rescaled densities at time t = are 

//(mi, . . . , m„ 0) = -i- f^'imi, m,), j = l,2,...,No (5.4) 

where the function /o has been defined in ()2.1ip . By ()5.4p . for every j > 1: 

lim ||//(0)-p^o/o'''llLi((M+y) = 0. (5.5) 

V, No^+oo 



In order to state our main result we first recall that, for all non negative initial data in 
the Cauchy problem for the coagulation equation (jl.ip with kernel equal to one has a unique non 
negative solution in C([0, +oo); L^(M+)) (cf.[I7], Theorem 2.1). 

Our main result is then the following: 

Theorem 5.1 Let {/j(t)}j^i be the solution of system 113. 13\) with initial data defined in 113. 5\ }. Then, 
{fj^ {'t)}f=i is the sequence of rescaled densities defined by /15.2\) : 

Vi>l: hm ||//(t)-/(t)®^||Li((M+).)=0 (5.6) 

Nr. 

where f is the unique solution in C([0, +oo); L^(M"'")) of the coagulation equation ( li. j|) with kernel 
K = 1 and initial datum pofo, /o given by \2.11\) . 

The proof of Theorem 15.11 is done in two steps. One is to prove that the sequence {//(i)}j^i 
converges to a sequence of functions {/j^(i)}jeN* that satisfy an infinite set of equations. The second 
is to prove that for all j > 1 and t > 0, fj^{t) = f{t)'^K The first uses the explicit expression of 
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the functions fj as a finite sum and a kind of dominated convergence. The second is done proving 
the uniqueness of solutions of the new infinite system of equations. We actually start by proving this 
uniqueness result. 

For the sake of notation let us introduce operators Wj defined as follows. Given a function 
(p G L^{{M.^y^^) we call Wj the operator such that Wj[(p] is: 

/•oo 

Wj[(p]inij) := VGj[ip]{nij) - j dnip{nij,fi), (5.7) 

Jo 

where Gj is defined by in (|2.14p . Therefore, using Lemma ()2.3p it is easily seen that Wj is a linear 
and continuous operator from L^((M"'")-'"''^) to L^{{'K~^y) whose norm satisfies: 

l|M^.II<f- (5.8) 



Lemma 5.2 Let f € C([0, +oo); L-'^(M"'')) be the unique non negative solution of the coagulation 
equation 

-j^ rm poo 

dtf{m,t) = - dfi fim-ii,t)f{i^,t)- fim,t) di^f{i^,t) (5.9) 



2 _ 

with initial datum /(0,m) = pofoi^^), fo given by 112. Then, the sequence of functions defined as 

/f(t,m,):=/(t)®^'(m,), j € N* (5.10) 
is the unique non negative solution of the system: 



i=i 



mi 



df°° 1 ^ 

-^{mj,t) = j dfi fj^-i^{mi,...,mi- fi,...,mj,ii,t) + 



poo 

3 dufjliimi,. . . ,mj,^i,t), forj = l,2,... (5.11) 
Jo 



with initial data f°°{0) = p^f^^ such that f°° G C([0, +cx)); Li(M+)J). For all t > 0, the sequence 
{/f (ObeN* satisfies: 

Vi>l: ||/r(i)llLi((R+).)<p'o- (5-12) 

Proof of Lemma 15.21 A straightforward calculation shows that the sequence {/?°(t)}jgN* defined in 
(|5.10p is indeed a solution of system (j5.1ip with initial data {p'f^f®-'{m.j)}j(z^*. 

On the other hand, since / G C([0, +cx)); L^(]R+)) satisfies we have / G C^((0, +00); L^(M+)) 
and then /°° G C([0, +00); Li(]R+y) n €^((0, +00); Li(M+)J). Due to the fact that {/f (t)}jeN* 
satisfies (|5.1ip . we obtain, after integration of the j— th equation over (M+)-': 

d j 
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and then, for all t > 0: 

\\fr(t)\\LHiR+y) < ii/r(o)iiiH(M+M) = Po- (5-13) 

That proves (l5J2]l . 

In order to prove uniqueness notice that the system (|5.1ip may be written using operator Wj as 



from where 



OOl 



/fW = /r(0) + fdhW, f°l,{t,) (5.14) 
Jo 



Using this same formula for fj^i, we obtain after M iterations: 



'3 
AI 



j.n 

= E ^^+1 • • • w^+n-i fru^) + 

n>0 

rt rtM 

+ dti... dtM+lWj . . .Wj+M fj^M+li^j+M+i;tM+l)- (5.15) 

Jo Jo 

Let us assume now that there exist two different non-negative solutions, {/ij,fe(i)}jgN* , A; = 1, 2, of 
([^TTD such that hj^k G C{[0, +oo); L^{R+y) with the same factorized initial datum {p^ (/o)®-'}jGN* • 
Then by ()5.15p and (|5.8p we deduce that for every M > 0: 

\\hj,l{t) - /ii,2(t)||Li{(R+y) < ^M + l)\ ^ 

x\\Wj . . . Wj+M {hj+M+l,l{tM+l) — ^j+M+l,2(iM+l)) llLi((R+)J+*-f+i) 

- (M + l)! V2 y (j _ 1)1 (,ll^+A^+l,l(*Af+lj||Li((R+).+M+i) + 

+ ll^i+Af+l,2(iAf+l)j|Li((R+)j+A/+i) j , (5.16) 



By our assumptions on {hj^k{'t)}ji=n*, k = 1,2, they both satisfy ()5.12p . Thus, plugging ()5.12p into 
(|5.16p we deduce: 



\h,,i{t) - /i,,2(t)||Li((R+).) < 2po (M + l)! U j (j - 1)! 



(5.17) 



Using ()3.1ip we obtain 

\\hj,i(t) - /i,,2(t)||Li((R+y) < (3^0*)"^+' (5.18) 
and then, for every t < l/{3po)- 

lLi((M+)i) 



h,,i{t) - hj,2{t)\\LHiR+y) < , lim {3pot)^'+' = (5.19) 
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and the uniqueness follows but only for t G [0, l/3po)- This argument may be now iterated as follows. 
Suppose that ()5.19p holds for t G [0, T) for some T > 0. Consider now (|5.16p for t G [T,T + l/3po) 
(since it holds for all t > 0). Since ()5.12p holds also true for all t > 0, it follows as before that (|5.19p 
also holds for t G [T, T + l/Spo ) • Therefore, global uniqueness is proven. This ends the proof of Lemma 
521 □ 



(5.20) 



Lemma 5.3 Under the same hypothesis than in Theorem \ 5. 1\ for allT > and all j > 1: 



lim sup ||//(t) - 9j{t)\\LmR+y) = 0. 

V,Afo-i'+ooo<t<r 



where 



g,{t, m,) = Pi'-'^I^W, o W,+, • • • o Ty,.+„_i[/f +^(m, 



n=0 



(5.21) 



is such that gj G C([0, +oo); Li(M+)J) n C^{{0,oo);L'^{R+y) and gj > 0. Moreover the sequence 
{9j}j€N* satisfies system 115.11]) for t G (0, oo). 

Proof of Lemma 15.31 Let us prove first that the series in (|5.2ip defines a function gj G 
C([0,Ti);Li(M+)J) n Ci((0,Ti);Li(]R+)J) for some n > 0. To this end we first deduce the follow- 
ing estimate from ()5.8p and the the explicit expression ()5.4p of the initial datum: 



Wj o ■ ■ ■ o 



Li{(M+)J) 
j{j + l)...{j + n-l). 



from where, using again (|3.1ip we deduce , for all j > 1, ?i > and t > 0: 



(5.22) 



j+n ^ 

Po 



WjO---0 Wj+n-\[fi 



Li((R+)J) 



<2^-V^o(3pot)". 



(5.23) 



It then follows that the series in ffOT]) defines a continuous function gj G C{[Q, ^);L^{{m+y)). A 
similar argument shows that gj G C^((0, ■^)] L} {{Ef^y)) and that 



dg,{t) 
dt 



4.11 — 1 



n=l 



By construction the sequence {fj {t)}^li satisfies, for every t > 0: 



//(i,m,) = e-^VrK;0)-/ dtie-^ 



d/i//+i(mj,/i;ti)+ (5.24) 



+ 



2 Jo 



dti e" 



E / d^i fjj^^{mi,. . . ,mi- IX,. . . ,m.j,iJ,;ti). 
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that, using the operator Wj defined by ()5.7p . can be rewritten as: 



No-j 



n>0 



(j+n-2)(j + n-l) 



(j + n— 1)/, , \ + ^ — I T7 

^ T^^.+„_,e- ' ' *"/]:;„(0), (5.25) 



where we set t_i = 0, to = ^• 

Property ()5.20p foUows from the two fohowing facts. The first is that the sum: 



n>0 





rtn-1 


/ dh.. 


■ / dtn 


Jo 


Jo 



Li({R+)J) 



is dominated by a convergent series uniformly for t € [0, l/6po). The second is that for every j > 1, 
every n > 1: 



V f f 1^ 3(3-1) ff , (j+n-2)(j + n-l) . ^ U + n-l)(j+n) ^ \ 

/0-5>+OO \j J 

= P'o^^^W^. ° • • • ° W.+n^AC^''] (5-26) 



y, Afo-!>+oo 



n! 

in Li((]R+)J), uniformly for t G [0,r] for any T > 0. 



We start proving the first. By (j5.4p and (|5.8p we obtain: 

o . . . o T^,+„_i[//+„(0)]||^,((j^^),) < (0 i(j + 1) . . . (j + n - l)||//+„(0)||ii((M+).) 



3\" /iyoy+" 



Since in the expression ()5.25p t„ < tn-i < ■ ■ ■ < < = \t follows that tk-i — > for any 
k = 1, . . . ,n and by (|5.27p we deduce: 



dtn 



Li((R+)J) 



j(j-l) ^-t ^ X (j + n-2)(j + n-l) X (j + n^l)(j + n) ^ 

. . . e 2v Wj+„_ie /j+„(0) 

J(j-l) ^ X (j+n-2)(j + »i-l) . X (j + n-l)(j + n) 

2V 1'' P 2V l,Cn-l l^njg 2V ' 



(2) + + (^-^j . (5.28) 
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Using again ()3.1ip and the hypothesis ()5.ip we deduce, for all j > 1, n > and t > 0: 

^'■l-^h f »,UTr (j+»-2)(j + n-l) . . (j + n-l)(j + n) , y 

Li((R+)i) 

<4.^-^{<opotY. (5.29) 

Estimate (j5.29p shows that the sum giving fj{t) in ()5.25p is dominated, uniformly for t G [0, ri), 
Ti = g^, by a convergent series. 

In order to deduce (j5.26p for every j > 1 and every n > 1 uniformly for t G [0, ri) we notice first of all 
that, as y ^ +oo: 

dine — 2v^(*-*i) • • • e (tn-i-t„)g — ^ / tft„ ee — (5.30) 

uniformly for t in any compact of [0,+oo). On the other hand by (|5.5p and the continuity of the 
operators Wk we have 



lim 

Nr. 



j+n-1 



/An(o) - pre-''' 



Li((R+y) 







(5.31) 



Prom (fOOD and (IOT]l . HOG]) follows and therefore: 



n>0 



in C([0,Ti);Li(M+)) n Ci((0, n); Li(]R+)). Since //(t) > for ah j G N* and t > it follows that 
gj(t) > for t G [0,ri) and all j. 

A straightforward computation shows that {gj}jeN* satisfies system ()5.1ip for t G [0,ri) and this 
proves the Lemma 15.31 for T G [0,ri). 

In order to extend this result for T > ri we first notice that for all f > 0: 



J(j-i) 



e-^V/(n/2) + / dh 



e 2v 



n>0 



dtr,e 2 V " 



{j + n-2)(j + n-l) , ^ {j + n-l)U + n) , 

^-it^-i-t.) e- ' 2v ^ *^//^Jri/2), (5.32) 

In order to pass to the limit as y ^ +cx3 and N^/V ^ in. (|5.32p we use the same two arguments 
as for ()5.25p . More precisely, using Lemma 221 we obtain: 

\\W, o . . . o ty,+„„i[//+Jri/2)]||^,((j^^^^.^ < j{3 + 1) . . . (i + n - l)||//+„(0)||^i((ij+y) 

< (0"i(j + 1) . . . (i + n - 1) ly^y^ . (5.33) 



. . . e 
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and then, arguing as before: 

(j + n-2)(j + n-l) ,^ . N (j + n-l)(j + 7i) , 



dtr 



e 



Li({I8+)J) 

< 4^"i(6/9ot)". (5.34) 



On the other hand, by the continuity of the operators Wk and the convergence (|5.20p that we have 
just proved for T € [0, ti) : 

nJ'^+oo ^^-+-1 [//+n(n/2) - 5,(ri/2)] = (5.35) 

From ()5.34p and ()5.35p we deduce that, as V ^ oo, Nq/V ^ po 

//(t + y) ^/i,(t) in C([0,ri);Li(IR+))nCi((0,ri);Li(M+)) 
= 5^ ;^/'o^"^^^ ° ■ ■ ■ ° ^i+-ife(n/2)] > Vt G [0, n). 

n>0 

We define now the function gj for t G [ri, 3ri/2) as : 

gj{t) = hj{t - ri/2), Vt e [ti,3ti/2). 
Since {5j(t)}jgAr* satisfies system ()5.1ip for t G [0, ri) it is such that: 

i + y) = E ;^'^o^"^:'- ° ■ ■ ■ ° ^:;+-ife(n/2)] yt G [0, n/2) 

n>0 

It then follows hj{t) = gj{t + ri/2) for all t G [0,ti/2). We deduce that gj G C([0, 3ti/2); Li(]R+)) n 
Ci((0,3ri/2);Li(M+)) and 

//(t)^5,(i) in C([0,3ri/2);Li(M+))nCi((0,3ri/2);Li(M+)) 

No 

as A/'o ^ oo, y ^> oo, -^^^ — > pq. 

Since, as a plain calculation shows again, the sequence {hj{t)}j^^* satisfies system ()5.1ip for t G [0,ti) 
and hj{0) = gj{Ti/2) the sequence {gj}^^^, satisfies system (jS.lip for t G [0,3ri/2). 

We have extended in that way the previous result on [0,ri) to the interval [0,3ri/2). This 
procedure may be repeated to obtain the result in all finite interval [0, T) from where Lemma 15.31 
follows. □ 



Proof of Theorem 15.11 From Lemma 5.3, for any T > 0, converges to gj in C([0, T); L^(]R^)-') 

and {gj}jeN* is a non negative solution of system ()5.1ip with initial data {/Og/^'^lieN* . By Lemma 
5.2 we have gj = for all j > 1 and Theorem 15.11 follows. □ 



24 



Remark 5.4 We have considered afactorized initial data for the finite system of particles (cf. i2.10\) ). 
We could also have considered a more general initial data of the form: 

P7Vo(m^o>0) = (^o)! *Wo(miVo), ^V(miv,0) = for N ^ Nq 

with 

^Noi^No) > 0, y drUNf, ^^Noi^No) = 1, (5-36) 
Then, for any j = 1, . . . , Nq, the correlation functions at time t = would have been 

/,(m„0) = f; I dm,+i...l dmNPNiniM,0) = (^^^ Fj''°\in,), (5.37) 

where Fj^°\inj) is the j -particle marginal distribution (in the sense of probability measures) of 
PNoiniNo), namely: 

Fj^°\mj) := j druj+i ... j dm.^o'^ NoimNo) 
and, by definition (see liS.clb]) ). we have: 

J dmj Fj^"\inj) = j drnj j dmj+i ... j dmNo"^ No{t^No) = 1- 



Then, the rescaled correlation functions at time t = would have satisfied 

//(0)||li{(k+)j) 



lim ||//(0)||ii((K+).)=/^o (5-38) 



Now, assuming 



for some fo{m) such that 



ivn 



lim ||F,^^"^(0)-/o®^||^i((M+).)=0 (5.39) 



fo{m) > 0, J dm /o(m) = 1, 



we could have proved exactly the same results that we got for the a priori factorized case. 
Choosing an intial datum as in 15. 39\) would have meant to assume propagation of chaos to hold at 
time t = 0. On the other hand, what we did has been to assume hypotheses of molecular chaos to hold 
at time t = 0. 
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